Landauer theory of ballistic torkances in non-collinear spin valves 
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We present a theory of voltage-induced spin-transfer torques in ballistic non-collinear spin valves. 
The torkance on one ferromagnetic layer is expressed in terms of scattering coefficients of the 
whole spin valve, in analogy to the Landauer conductance formula. The theory is applied to 
Co/Cu/Ni(001)-based systems where long-range oscillations of the Ni-torkance as a function of 
Ni thickness are predicted. The oscillations represent a novel quantum size effect due to the 
non-collinear magnetic structure. The oscillatory behavior of the torkance contrasts a thickness- 
independent trend of the conductance. 

PACS numbers: 72.25.Mk, 72.25.Pn, 75.60.Jk, 85.75.-d 



I. INTRODUCTION 

The predictionii^ and realization^ of current-induced 
switching of magnetization direction in epitaxial mag- 
netic multilayers stimulated huge research activity re- 
lated to high-density writing of information. The 
simplest systems for this purpose are spin valves 
NM/FM1/NM/FM2/NM with two ferromagnetic (FM) 
layers (FMl, FM2) separated by a non-magnetic (NM) 
spacer layer and attached to scmiinfinite NM metallic 
leads. The electric current perpendicular to the layers 
becomes spin-polarized on passing the FMl layer with 
a fixed magnetization direction. In non-collinear spin 
valves, subsequent reflection and transmission of spin- 
polarized electrons at the FM2 layer results in a spin 
torque acting on its magnetization the direction of which 
can thus be changed. Majority of existing experimental 
and theoretical studies of these spin-transfer torques re- 
fer to a diffusive regime of electron transport in metallic 
systems, see Ref. H for a review. 

Magnetic tunnel junctions with the NM spacer layer 
replaced by an insulating barrier have attracted atten- 
tion only very recently; in these systems voltage-driven 
spin-transfer torques^ as well as effects of finite bias^i^ 
can be studied. The concept of torkance. defined in the 
small-bias limit as a ratio of the spin-transfer torque and 
the applied voltage^ represents an analogy to the con- 
ductance. It becomes important also for all-metallic spin 
valves with ultrathin layers2ii£ where a ballistic regime 
of electron transport can be reaHzed. 

The latter regime is amenable to fully microscopic, 
quantum-mechanical treatments. All existing theoretical 
approaches to the torkance, both on model^i^ '^^i^^ and ab 
initio^'^ levels, are based on a linear response of various 
local quantities inside the spin valve to the applied bias. 
The local quantities used range from scattering coeffi- 
cients of the individual layersi^ over local spin currents^ 
to site- and orbital-resolved elements of a onc-particle 
density matrix.— These methods contrast the well-known 
Landauer picture of the ballistic conductanc o^^i^^ which 
employs only transmission coefficients between propagat- 



ing states of the two leads. 

In this paper, we present an alternative theoretical ap- 
proach to ballistic torkances that yields a result simi- 
lar to the Landauer conductance formula, i.e., we relate 
the torkance to scattering coefficients of the whole spin 
valve. This unified theory of both transport quantities is 
used to discuss a special consequence of ballistic trans- 
port, namely, a predicted oscillatory dependence on Ni 
thickness in a Cu/Co/Cu/Ni/Cu(001) system. The pre- 
sented study reveals a relation between the torkance and 
the properties of individual parts of the spin valve which 
might be relevant for design of new systems. 



II. THEORY 

A. Model of the spin valve 

Our approach is based on an effective one-electron 
Hamiltonian of the NM/FM1/NM/FM2/NM system. 



H = Hq + 7ini • (T + 72n2 ■ cr. 



(1) 



where Hq comprises all spin-independent terms, 71 = 
71 (r) and 72 = 72(1") denote magnitudes of exchange 
splittings of the FMl and FM2 layers, respectively, ni 
and n2 are unit vectors parallel to directions of the ex- 
change splittings, and the cr = {crx,cry,(Jz) are the Pauli 
spin matrices. The angle between ni and n.2 is denoted 
as 9. The spin torque r is defined as time derivative of 
the electron spin, represented (in units of Bohr magne- 
ton) by operator cr. This yields (with h = 1) the total 
spin torque as 



r = -i [cr, H] = T1+T2, 
where the quantities 

T,=27jnjXcr, J = 1,2, 



(2) 



(3) 



can be interpreted as torques experienced by the two FM 
layers. Obviously, the torque Tj is perpendicular to the 
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vector rij and it can thus be decomposed with respect to 
the common plane of the two magnetization vectors into 
its in-plane (t,||) and out-of-plane (Tj_l) components, see 
Fig. 1 for j = 2. The unit normal vector of the plane is 
given by 1/ = rii x 112/ sin 6. 




FIG. 1: (Color online) The in-plane (t2||) and out-of-plane 
(t2x) components of the torque T2 = T2\\ + T2± experienced 
by the FM2 layer. For details, see text. 



B. In-plane torkance 



The basic idea for the in-plane torkance on the FM2 
layer rests on the orthogonality relations rij ■ Tj = 0, 
j = 1,2, from which the size of T2\\ can be written as 



(n2 X u) ■T2 



ni ■ T2 



ni 



(4) 



see Fig. 1. The total torque r, being a full time derivative 
of cr, in (HI) plays a key role in the following treatment. 

Our approach applies to systems consisting of the left 
{£) and the right {TZ) semiinfinite NM leads with an 
intermediate region (X) in between; the latter contains 
both FM layers and the NM spacer of the spin valve. 
Projection operators on these regions are denoted re- 
spectively as n^, Il-ji and Tlx; they are mutually or- 
thogonal and satisfy He + Hj + H-jz = 1. The Hamil- 
tonian ([T]) is assumed to be short-ranged (tight-binding), 
not coupling the two leads, i.e., HcHH-jz = 0. The 
leads are in thermodynamic equilibrium at zero temper- 
ature. A general linear response theory can be formu- 
lated for a Hermitean operator Q = that is local, 
i.e., not coupling neighboring parts of the system, so that 
Q = UcQUc + rtiQIIi + n-jiQ^n- Its time derivative 



D = -i[Q,H] 



(5) 



is assumed to be localized in T, i.e., D = HxDYlj. These 
properties make it possible to remove the semiinfinite 
leads from the formalism. 

The resulting response coefficient describing the 
change SD of the thermodynamic average of the quantity 
D due to an infinitesimal variation Sfic of the chemical 
potential (Fermi energy) of the C lead is given by 



SD_ 

Sue 



1 

2^ 



Tr {Q (Ft^GT^G" - F^G'T^G'^)} , (6) 



where the trace (Tr) and all symbols on the r.h.s. are 
defined on the Hilbert space of the intermediate region X, 
in particular the Q in ([6]) abbreviates HjQUx- The other 
symbols in ^ refer to the antihermitean part of the C 
and TZ selfenergies, Tc,n{E) = i [S^ ^^(E') 
and to the retarded and advanced propagators 



(7) 



where Y.''^''{E) = S^^l-^) + denotes the total 

selfenergy. Omitted energy arguments in ^ are equal to 
the Fermi energy of the equilibrium system [E = Ep). 

The proof of ([6]) is based on non-equilibrium Green's 
functions (NGF) for stationary statesi^ and it is similar 
to a previous derivation in Ref. [igI . The starting point is 
an expression for the variation of Z), 



1 r°° 

6D = — Tr{G''{E)DG''{E)6J:<{E)}dE, 

J — oc 



(8) 



where the variation of the lesser part of the selfenergy at 
zero temperature is given by 



S^<iE) ^ d{E - EF)Tc{E)S^ic. 



(9) 



The assumed properties of H, Q and D lead to a com- 
mutation rule for the selfenergy. 



(10) 



which is proved in the Appendix and which in turn yields 
a relation 



G''{E)DG''[E) 



i[G''{E)Q~QG^{E)] 
+ G"'[E)QT[E)G'-{E), (11) 



where T{E)=i [17 {E) - {E)] = Tc (E) + rn{E). The 
result ([6]) follows then from an identity for the spectral 
density operator, 

i [G''{E) - G^iE)] ^ G''{E)T{E)G'-{E). (12) 

Note that the final response coefficient ((6)) obeys a perfect 
C-TZ symmetry, i.e., SD/S^i-ji = —SD/S^c- 

It should be emphasized that the derived general re- 
sult dS]) and its perfect C~TZ symmetry are valid only for 
operators D that can be formulated as a time derivative 
of a local operator Q according to ((S]). In the present 
context of spin valves, this is the case of the usual parti- 
cle conductance and of the in-plane torkance (see below). 
The out-of-plane torkance requires a different approach 
based on the more general relation ([8|) , see Section III C| 
its symmetry properties for symmetric spin valves were 
discussed in details, e.g., in Ref. [TtI. 

Application of the derived formula ([6]) to the transport 
properties of the spin valves is now straightforward. The 
usual particle conductance G is based on the operator Q 
being a projector on a half-space containing, e.g., the TZ 
lead and an adjacent part of the X region. This results 
in the well-known expressior>i^ 



G 



2tt 



Tr (FtjGT^G'^) , 



(13) 
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where atomic units (e = ?i = 1) are used. The in-plane 
torkance C|| on FM2 according to ^ is obtained from 
Q = rii • cr in ([6]). This yields C\\ = Ci/ sin0, where 



Ci 



1 

2^ 



Tr {m • a- {VnG^TcG-' - r^GT^G")} . (14) 



The two terms on r.h.s. can be related to spin fluxes on 
two sides of the FM2 layer. The expression rep- 
resents our central result. The operators r£ and V-ji 
are localized in narrow regions at the interfaces C/I and 
respectively. The Green's functions (propagators) 
for points deep inside the spin valve thus enter neither 
the conductance ((T3)) . nor the in-plane torkance (fT4)) . 



C. Out-of-plane torkance 

A similar approach for the out-of-plane torkance on the 
FM2 layer employs an infinitesimal variation 5xy2 of its 
magnetization direction due to a variation 59 of the angle. 
The FMl magnetization direction as well as the plane of 
the two directions ni, n2 remain fixed, i.e., 5ni ~ 5v 
This leads to 5r\.2 ^ v y. r\.256 and from ([1]) also to 



H' 



6H 



1 



72(1^ X n2) ■ (T ^ -V> ■ T2, 



(15) 



so that the size of T2± coincides (up to factor of 2) with 
angular derivative of the effective Hamiltonian H. The 
NGF formulation of the out-of-plane torkance rests on re- 
lation © applied to the operator D = 2H' , see (fT5|) . with 
variation of the selfenergy STi^(E) due to an infinitesimal 
variation of the chemical potential (5/i£ given by ^ and 
similarly for (5S<(iJ) due to the Sfi-ji. This yields then 
response coefficients Cc = SD/Sfj-c and C-ji = SD/Sfi-jz 
for the out-of-plane torque with respect to chemical po- 
tentials of the C and TZ leads expressed as 



1 



G£,K = -Tr(i/CT£,KG''), 



(16) 



By employing a simple consequence of p^ . G° = (1-1- 
iG"r)G'', cyclic invariance of trace, angular indepen- 
dence of the selfenergy of NM leads, T,"' = E'" = 0, and 
the rule G^H'G^ ~ G"", the response coefficients can be 
recast into 

Cc.n = -Tr {G"Tc.n [1 + iG^iVc + F^)]} , (17) 

which contain again only propagators at points close to 
the C/T and T/TZ interfaces, similarly to (fT51 [Ti|) . 

The applied bias has to be identified with the difference 
fJ'C ~ fJ'K and the out-of-plane torkance on FM2 is thus 
given by C± = {Cc — C-ji)/2. Since the Hamiltonian ([T|) 
does not contain spin-orbit interaction, the spin reference 
system can be chosen such that both unit vectors ni and 
n2 lie in the x — z plane. This implies that H is essentially 
time-inversion invariant and it can be represented by a 
symmetric matrix, — H\ the related quantities G*"'" 



and T cTi are symmetric as well. As a consequence, the 
transmission- like terms in Cc and C-jz are the same, i.e., 
Tr (G'T^C^Ftj) = Tr(G"'^7^G'^^£). The resulting out- 
of-plane torkance 

Cj. - ^Tr {G'^[Tc{l + iG-^F^) - Fk(1 + iC^Tn)]} 

(18) 

contains thus only reflection-like terms. 



D. Landauer formalism 

The Green's function expression for the conductance 
(fT^ can be translated in the language of scattering 
theory^ the counterparts of the torkances p4l [T8)) are 
interesting as well. In the present case, propagating 
states in the C and TZ lead will be labelled by A and 
p, respectively. Moreover, a spin index s =t, J, has to 
be used even for NM leads, since non-coUinearity of the 
spin valve gives rise to full spin dependence of scattering 
coefficients. 

The conductance (|13p is given by the Landauer formula 

G = (Stt)"^ X^Apss' , where <ps',As denotes the 
transmission coefficient from an incoming state As into 
an outgoing state ps' The in-plane torkance coefficient 
(fT4l) can be written as 



^1 ^ 27r ^ {ni-(T)s"s' {tps',\st*ps'^^Xs - t\s',pst*Xs",ps) > 
Xpss' s" 

(19) 

whereas the out-of-plane torkance ((T8)) can be trans- 
formed into 



X' s'Xs X' s' .Xs 



5^ ^ P's'.ps^p' 



s' ,ps I ' 



, XsX' 



(20) 

where rA'^'.As {rp's',ps) denote reflection coefficients be- 
tween states of the C (TZ) lead. This result represents 
analogy to the Landauer formula and it completes the 
unified theory of conductances and torkances. 



III. RESULTS FOR CU/CO/CU/NI/CU(001) 
AND THEIR DISCUSSION 

The developed formalism allows to study properties of 
spin valves with ultrathin layers, which is yet an experi- 
mentally unexplored area; here we demonstrate its use for 
understanding unexpected features of ab initio results. 
The results discussed below were obtained using the re- 
sponse of spin currents on both sides of the FM2 layer 
implemented within the scalar relativistic tight-binding 
linear muffin-tin orbital (TB-LMTO) methodi^i^ sim- 
ilarly to our previous transport studiesJ^i^i As a case 
study, spin valves Cu/Co/Cu/Ni/Cu(001) with face- 
centered cubic (fee) structure were chosen. All atomic 
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positions were given by an ideal fee Co lattice with sharp 
interfaces between the neighboring FM and NM layers. 
The spin valves discussed below consist of a Co layer of 
5 monolayer (ML) thickness separated by a 10 ML thick 
Cu spacer from a Ni layer of varying thickness, embed- 
ded between two semiinfinite Cu leads. Self-consistent 
calculations within the local spin-density approximation 
(LSDA) were performed only for coUinear spin valves 
{6 = or 6 = tt) while the electronic structure of 
non-coUinear systems was obtained by rotation of the 
exchange-split potentials of the Co and Ni FM layers. 
Particular attention has been paid to the convergence 
of torkances with respect to the number of k|| vectors 
sampling the two-dimensional Brillouin zone (BZ) of the 
system; in agreement with Ref. US we found that reliable 
values of the out-of-plane torkances require finer meshes 
than for the in-plane torkances. The presented data were 
obtained with 6400 kn points in the whole BZ. 
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FIG. 2: (Color online) Calculated transport coefficients (per 
interface atom) as functions of Ni thickness: (a) the conduc- 
tance (C) for three values of the angle 6 and the in-plane 
(C||) and out-of-plane (Cx) Ni-torkances for 6 = tt/2 in spin 
valves Cu/Co/Cu/Ni/Cu(001), (b) the real (Re) and imag- 
inary (Im) parts of the spin-mixing conductance (C""") of 
Cu/Ni/Cu(001) systems. 

Figure 2a displays the calculated conductances for par- 
allel {9 = 0); antiparallel {9 = tt) and perpendicular 
{6 — tt/2) orientations as well as Ni-torkanccs in the 
latter case as functions of Ni thickness. The most pro- 
nounced feature of the transport coefficients are oscilla- 
tions with a period of about 12 ML seen in both com- 
ponents of the torkance. These oscillations reflect the 
perfect ballistic regime of electron transport across the 
whole spin valve. In addition, they contradict a generally 
accepted idea of very short magnetic coherence lengths of 



a few interatomic spacings, or, equivalently, of the spin- 
transfer torques as an interface propertyi^i^ i '^'^ Very re- 
cently, spin-transfer torques in antiferromagnetic metal- 
lic FeMn layers have been investigated theoretically 
it has been shown that the torques are not localized to 
the interface but are effective over the whole FeMn layer. 
However, no oscillatory behavior of the total torkance as 
a function of the layer thickness has been reported. The 
nature of the predicted oscillations deserves thus detailed 
analysis, including also a discussion of their stability with 
respect to structural imperfections and of their absence 
in the conductance (see Fig. 2a). 

Oscillations similar to those in Fig. 2a have recently 
been obtained for a different quantity of a simpler sys- 
tem, namely for the spin-mixing conductance C"""^ of 
epitaxial fee (001) Ni thin films attached to Cu leadsJ^ 
The real and imaginary parts of the complex C'™"' are 
related to two components of the spin torque experi- 
enced by the FM film due to a spin accumulation in one 
of the NM leadsr the calculated values of C^™^ for the 
Cu/Ni/Cu(001) system arc shown in Fig. 2b. The oscil- 
lation periods of the torkance and the spin-mixing con- 
ductance are identical which indicates a common origin 
of both. The physical mechanism behind the C™"^ oscil- 
lations was identified with an interference effect between 
spin-f electrons propagating across the Ni film from the 
TZ lead to the £ lead and spin-| electrons propagating 
backwards. This effect is expressed by a spin-mixing term 
- tr{rTzG\TcGl) in the C™'', where the Gl'" {s =t i) 
denote spin-resolved propagators and the trace (tr) does 
not involve the spin index.— The particular value of the 
oscillation period follows from a special shape of the spin- 
polarized Fermi surface of bulk fee NiM- 

The oscillations of C*™"^ have been found fairly sta- 
ble against Cu-Ni interdiffusion at the interfaces'^ the 
same stability can be thus expected for the torkance os- 
cillations in the spin valve. The relative stability can 
be understood as an effect of the large oscillation period 
(~ 12 ML): intermixing confined to a very few atomic 
planes at interfaces reduces the oscillation amplitude 
rather weakly. This feature contrasts, e.g., sensitivity 
of the interlayer exchange coupling in magnetic multilay- 
ers mediated by a NM Cu(OOl) spacer with oscillation 
periods of ~ 2.5 ML and 6 ML, where even a very small 
amount of interface disorder reduces strongly especially 
the amplitude of the short period oscillationsj^ 

Another important aspect of the oscillations of the 
spin-transfer torques concerns their dependence on the 
thickness of the polarizing Co layer since ultrathin lay- 
ers in general might amplify ballistic and interference ef- 
fects. Figure 3 presents the Ni-torkances in the same 
Cu/Co/Cu/Ni/Cu(001) spin valves calculated for three 
different Co thicknesses, namely 5, 15 and 25 ML. It can 
be seen that the oscillations persist and have the same pe- 
riod in all three cases. Their amplitudes depend slightly 
on the Co thickness: the initial increase from 5 to 15 Co 
ML is accompanied by a small reduction (of about 20 
%) of the amplitudes whereas further increase from 15 
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FIG. 3: (Color online) Calculated in-plane (Cy) and out-of- 
plane (C_l) Ni-torkances (per interface atom) as functions of 
Ni thickness in spin valves Cu/Co/Cu/Ni/Cu(001) for 6 = 
7r/2, 10 ML Cu spacer and for three different Co thicknesses. 



IV. CONCLUSION 

We have addressed two important aspects of spin- 
transfer torques in non-coUinear spin valves with ultra- 
thin layers. First, we have shown that the in-plane and 
out-of-planc torkance on one FM layer can be expressed 
by means of the transmission and reflection coefficients, 
respectively, of the whole spin valve, in close analogy to 
the Landauer formula for the ballistic conductance. Sec- 
ond, a novel oscillatory behavior for Ni-based systems 
has been predicted due to the mixed spin channels. The 
oscillations with Ni thickness are reasonably stable with 
respect to interface imperfections of real samples; how- 
ever, they are not present in the conductance but can 
be observed only in the Ni-torkance. The torkance os- 
cillations prove that the spin-transfer torques in ballistic 
spin valves are closely related to properties of their com- 
ponents, in particular to the spin- mixing conductances of 
individual ferromagnetic layers. 



to 25 Co ML docs not influence them appreciably. More 
detailed investigation of the effect of the thickness of the 
polarizing Co layer, including also the limiting case of 
spin valves FM1/NM/FM2/NM with a semiinfinitc po- 
larizing FM lead,— is beyond the scope of the present 
study. 

Let us now discuss the absence of oscillations in the 
conductance (Fig. 2a). We introduce propagators Gj" 
of an auxihary system NM/NM1/NM/FM2/NM, where 
NMl denotes the FMl layer with null exchange splitting. 
These propagators satisfy G"''" = + Gj °T'-[^'"G2 " 
where the denotes the i-matrix corresponding to the 
FMl exchange splitting 71 in ([1]). The conductance (fTS]) 
can be then rewritten as G = (27r)-iTr(r7?,G^A£G^), 
where Ac = {I + G5)r£(l + GfTf ) represents an op- 
erator localized at the FMl layer, i.e., at the left end of 
the FM2 layer. The latter trace can be most easily evalu- 
ated using the spin quantization axis parallel to the FM2 
magnetization direction n2. Since the propagators Gj" 
are now diagonal in the spin index s and the operator T-ji 
is spin- independent, the conductance does not contain 
spin-mixing terms, i.e., terms ~ tr(r-R,G2 ^ A^^ss/Gj ^z) 
for s ^ s' that result in interference effects involving 
different spin channels. For the torkance, however, the 
extra factor ni • cr in (fn|) provides the necessary spin 
mixing responsible for the oscillations, in full analogy to 
oscillations of the spin-mixing conductance. 

A recent study of spin-transfer torques in a tunnel 
jrmction Cu/Fe/MgO/Fe/Cu has predicted torkance and 
conductance oscillations with Fe thickness with a period 
^ 2 ML 4 These oscillations were ascribed to quantum 
well states in the majority spin of the Fe layer, i.e., to 
interference effects in a single spin channel. The oscilla- 
tions in the Cu/Co/Cu/Ni/Cu system — manifested only 
in the torkance — have thus a clearly different origin. 
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APPENDIX: PROOF OF THE COMMUTATION 
RULE FOR SELFENERGY 

The proof of the commutation rule (fTU)) rests on as- 
sumed properties of the operators H, Q, D (deflned in 
the Hilbcrt space of the total system) and of the projec- 
tors n£, Ht^, Hi, see the beginning of Section lll Bl Let as 
abbreviate projections of any operator X {X = H,Q,D) 
as HxXHj = Xjx, HxXHc = Xjc, etc. The assumed 
property of Q, namely Q = Qcc + Qix + Qmz, a con- 
sequence of the localization of D in X, namely Dxc = 
and Dci = 0, and the orthogonality of the projectors 
n^, Ht^, Hx lead to identities 

QiiHxc ~ HjcQcc, QccHci = HciQxx- (A.l) 

Similarly, a commutation rule 

QccHcc ~ HccQcc (A. 2) 

can easily be obtained from Dec = 0. 
The left selfenergy is given explicitly by 

Y.Y{E)^HxcQTiE)Hci, (A.3) 

where the Q2°'{E) denotes the retarded and advanced 
propagator of the isolated left lead. The relation (|A.2p 
implies immediately a commutation rule 

QccGTiE) = gTiE)Qcc; (A.4) 
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its application together with (jA.lilA.3|) leads to identities 

= HxcQccGc°'iE)Hcx = HicGc°'iE)QccHcx 
= HTcGTiE)HcxQxx^^TiE)Qix, (A.5) 



which are equivalent to the commutation rule (jlOp for 
the left selfenergy. The proof for the right sclfencrgy is 
similar and therefore omitted. 



* Electronic address: 'carva@karlov . mff . cuni.cz' 

t Electronic address: turek@ipm.cz. Also at Institute of 
Physics of Materials, Academy of Sciences of the Czech 
Republic, CZ-61662 Brno, Czech Republic 

^ J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996). 

^ L. Berger, Phys. Rev. B 54, 9353 (1996). 

^ E. B. Myers, D. C. Ralph, J. A. Katine, R. N. Louie, and 
R. A. Buhrman, Science 285, 867 (1999). 

* J. A. Katine, F. J. Albert, R. A. Buhrman, E. B. Myers, 
and D. C. Ralph, Phys. Rev. Lett. 84, 3149 (2000). 

^ A. Brataas, G. E. W. Bauer, and P. J. Kelly, Phys. Rep. 

427, 157 (2006). 
® J. C. Slonczewski and J. Z. Sun, J. Magn. Magn. Mater. 

310, 169 (2007). 

1. Theodonis, N. Kioussis, A. Kalitsov, M. Chshiev, and 
W. H. Butler, Phys. Rev. Lett. 97, 237205 (2006). 
® C. Heiliger and M. D. Stiles, Phys. Rev. Lett. 100, 186805 
(2008). 

^ D. M. Edwards, F. Federici, J. Mathon, and A. Umerski, 

Phys. Rev. B 71, 054407 (2005). 
° S. Wang, Y. Xu, and K. Xia, Phys. Rev. B 77, 184430 

(2008). 

^ P. M. Haney, D. Waldron, R. A. Duine, A. S. Nunez, 
H. Quo, and A. H. MacDonald, Phys. Rev. B 76, 024404 
(2007). 

^ X. Waintal, E. B. Myers, P. W. Brouwer, and D. C. Ralph, 

Phys. Rev. B 62, 12317 (2000). 
^ R. Landauer, Philosophical Magazine 21, 863 (1970). 



S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge University Press, 1995). 

H. Hang and A.-P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors (Springer, Berlin, 1996). 
K. Carva and I. Turek, Phys. Rev. B 76, 104409 (2007). 
D. C. Ralph and M. D. Stiles, J. Magn. Magn. Mater. 320, 
1190 (2008). 

L. Arrachea and M. Moskalets, Phys. Rev. B 74, 245322 
(2006). 

O. K. Andersen and O. Jepsen, Phys. Rev. Lett. 53, 2571 
(1984). 

^° L Turek, V. Drchal, J. Kudrnovsky, M. Sob, and P. Wein- 
berger, Electronic Structure of Disordered Alloys, Surfaces 
and Interfaces (Kluwer, Boston, 1997). 
K. Carva, I. Turek, J. Kudrnovsky, and O. Bengone, Phys. 
Rev. B 73, 144421 (2006). 

C. Heiliger, M. Czerner, B. Y. Yavorsky, L Mertig, and 

M. D. Stiles, J. Appl. Phys. 103, 07A709 (2008). 

M. D. Stiles and A. Zangwill, Phys. Rev. B 66, 014407 

(2002). 

Y. Xu, S. Wang, and K. Xia, Phys. Rev. Lett. 100, 226602 
(2008). 

^■■^ K. Carva and I. Turek, Physica Status Solidi A 205, 1805 
(2008). 

J. Kudrnovsky, V. Drchal, I. Turek, M. Sob, and P. Wein- 
berger, Phys. Rev. B 53, 5125 (1996). 



